home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CPODI.z / CPODI
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCPPPPOOOODDDDIIII((((3333FFFF))))                                                            CCCCPPPPOOOODDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CPODI   - CPODI computes the determinant and inverse of a certain complex
  10.      Hermitian positive definite matrix (see below) using the factors computed
  11.      by CPOCO, CPOFA or CQRDC.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE CPODI(A,LDA,N,DET,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAA COMPLEX(LDA, N)
  21.         the output  A  from CPOCO or CPOFA
  22.         or the output  X  from CQRDC.
  23.  
  24.      LLLLDDDDAAAA INTEGER
  25.         the leading dimension of the array  A .
  26.  
  27.      NNNN INTEGER
  28.         the order of the matrix  A .
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         = 11   both determinant and inverse.
  32.         = 01   inverse only.
  33.         = 10   determinant only.  On Return
  34.  
  35.      AAAA If CPOCO or CPOFA was used to factor  A  then
  36.         CPODI produces the upper half of INVERSE(A) .
  37.         If CQRDC was used to decompose  X  then
  38.         CPODI produces the upper half of INVERSE(CTRANS(X)*X)
  39.         where CTRANS(X) is the conjugate transpose.
  40.         Elements of  A  below the diagonal are unchanged.
  41.         If the units digit of JOB is zero,  A  is unchanged.
  42.  
  43.      DDDDEEEETTTT REAL(2)
  44.         determinant of  A  or of  CTRANS(X)*X  if requested.
  45.         Otherwise not referenced.
  46.         Determinant = DET(1) * 10.0**DET(2)
  47.         with  1.0 .LE. DET(1) .LT. 10.0
  48.         or  DET(1) .EQ. 0.0 .  Error Condition a division by zero will occur
  49.      if the input factor contains a zero on the diagonal and the inverse is
  50.      requested.  It will not occur if the subroutines are called correctly and
  51.      if CPOCO or CPOFA has set INFO .EQ. 0 .  LINPACK.  This version dated
  52.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  53.      Subroutines and Functions BLAS CAXPY,CSCAL Fortran CONJG,MOD,REAL
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.